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Abstract 

, We use a large-scale QTiQ calculation for B Li with microscopically derived two- 

body interaction to construct the OhQ Op-shell effective hamiltonian, electric 
quadrupole, and magnetic dipole operators. While the E2 and Ml QhQ opera- 
tors are one-body operators with free nucleon charges, the effective operators 
Q . are two-body operators with substantially different renormalization for the 

lf~) \ isoscalar and isovector matrix elements, especially for the E2 operator. We 

show that these operators can be very well approximated by one-body oper- 
J> ■ ators provided that effective proton and neutron charges are used. The ob- 

tained effective charges are compatible with those used in phenomenological 
^ shell-model studies. The two-body part of the effective operators is estimated. 



X 



Considerable effort has been devoted to derive the effective interaction used in the shell- 
model calculations from the nucleon-nucleon interaction [0-0]. On the other hand, much 
less work has been done to understand the effective operators employed in calculating differ- 
ent nuclear, usually electromagnetic, properties. In particular, a microscopic derivation of 
effective operators has been only partially succesful. It is well known that effective proton 
and neutron charges must be employed to describe the E2 transitions and moments. These 
charges are quite different from the free nucleon charges, typically the values of ~ 1.5e 
and e° ff ~ 0.5e are obtained for both light and heavy nuclei. Attempts to derive these 
charges microscopically, usually by perturbation ||, or by an "expanded shell-model" ap- 
proach H], yielded much smaller values. It should be noted that these effective charges 
correspond to a severely truncated single-major-shell, or OhQ, space. The question arises as 
to what causes the nucleon properties to change so significantly, is it mostly the result of 
the space truncation or the fact that nucleons are bound? Also the non-nucleonic degrees 
of freedom may play an important role. Other interesting questions are: how important are 
the higher than one-body parts of the effective operators and what is the j-dependence of 
effective charges? 

In this contribution we investigate how severe space truncation affects the electromag- 
netic operators. We use a large-space 6hQ shell-model calculation for 6 Li, with a micro- 
scopically derived two-body interaction, to construct an effective hamiltonian and effective 
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electromagnetic operators, which will exactly reproduce the 6h£l results in the Op-shell for 
the (0s) 4 (0p) 2 dominated states. This enables us to compare the resulting effective opera- 
tors with the bare one-body Op-shell operators and to extract the relevant effective charges, 
which allow us to determine the amount of renormalization, to study their j-dependence, 
and, eventually, to quantify the two-body content of the effective operators. Also we per- 
form the same derivation from the corresponding AhQ calculation to study the dependence 
on the space size and compare the rate of convergence for the effective hamiltonian and the 
effective operators. 

Recently, large-basis no-core shell- model calculations have been performed [^],|5| . In these 
calculations all nucleons are active, which simplifies the effective interaction, as no hole states 
are present. In the approach taken, the effective interaction is determined microscopically 
from the nucleon-nucleon interaction for a system of two nucleons and subsequently used 
in the many-particle calculations. To take into account a part of the many-body effects, a 
multi- valued effective interaction approach was introduced ||, which uses different sets of 
the effective interaction for different hfl excitations. In the latest application of the no-core 
approach, we derived starting-energy-independent hermitian two-body effective interactions 
from the Reid 93 nucleon-nucleon potential and applied them in the multi-valued ap- 
proach to A = 3 — 6 nuclei ||. In this study we use the results of this calculation for 6 Li 
presented in the third column of Table IV of Ref. ||. The many-particle calculation was 
done using the Many-Fermion-Dynamics Shell- Model Code |J in the m-scheme with dimen- 
sions approaching 2 x 10 5 . As in the previous large-scale no-core shell- model calculations 
|5|||, a reasonable description of the electromagnetic properties has been achieved using free 
nucleon charges. Our aim here is to study the renormalization of these operators, when the 
model space is severely truncated. 

For the OhQ dominated states of 6 Li shown in Table IV of Ref. |J it is possible to for- 
mulate an equivalent description purely in the Op-shell. We may divide the basis states of 
the 6hQ calculation into two subspaces, using the projectors P and Q, P + Q — 1. Here 
the P-space is spanned by the states |(0s) 4 (0p) 2 ). There are 10 such states in the Mj = 
m-scheme calculation and 8 in the Mj = 1 calculation, respectively. The Q-space is then 
formed by the rest of the almost 200,000 states. The effective P-space hamiltonian may 
be constructed by employing the Lee-Suzuki starting-energy independent similarity trans- 
formation method |TI| , which gives the effective hamiltonian PH e ^P = PHP + PHQuP, 
with the transformation operator to satisfying u = QuoP. In the case when the full space 
eigenvectors are known, like in our situation, this operator may be obtained directly. Let us 
denote the P-space states as \ap), and those which belong to the Q-space, as |«q). Then the 
Q-space components of the eigenvector \k) of the full-space hamiltonian can be expressed as 
a combination of the P-space components with the help of the operator uj 

(a Q \k) = J2(a Q \uj\ap} (a P | k) . (1) 

OCp 

If the dimension of the model space is dp, we may choose a set K of dp eigenevectors \k), for 
which the relation ([!]) will be satisfied. In our case we choose those states, which have the 
dominant OhQ component. Under the condition that the dp x dp matrix (atp\k) for \k) e /C 
is invertible, which is satisfied in the present application, the operator uo can be determined 
from (|l|). Consequently, the effective hamiltonian can be constructed as follows 
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{'j P \H eS \a P ) = [(lp\k)E k (k\a P ) 

+ J2( r rp\k}E k (k\a Q )(a Q \uj\ap) 



(2) 



It should be noted that P\k) = J2a P \oip)(a P \k) for \k) G JC is a right eigenvector of @ with 
the eigenvalue E k . The hamiltonian (0) is, in general, non- hermitian, or more accurately 
quasi-hermitian. It can be hermitized by a similarity transformation, which is determined 
from the metric operator P(l + u'u)P. The hermitian hamiltonian is then given by |TT] 



H. 



cff 



H, 



off 



(3) 



Similarly, a corresponding effective operator O c s can be constructed for any full space, 
e.g., electromagnetic, operator O so that it exactly reproduces the full space results for the 
P-space eigenstates. A double-similarity transformation |TJ| leads to the P-space operator 

(P + P^Q)6(P + 
(131) is then obtained 



associated with the hamiltonian (0) in the form O e s = P(l + uj^ui)P 
QujP). The operator associated with the hermitian P-space hamiltonian 
as EM31 



O 



off 



P(l + a; t a;)Pl 2 (P + P^Q)6(P + QwP) 



x 



P(l + cA;)P 



(4) 



Using the Eqs. (|1],|2|,|3|) , we constructed the effective hamiltonian, whose matrix elements, 
after performing the transformation from m-scheme to JT basis, are presented in Table |. 
In the same Table the well-known Cohen-Kurath matrix elements are shown. These 
were obtained by a least-square fit to experimental binding energies, relative to 4 He, and 
excitation energies of A = 6 — 16 nuclei. To make a meaningful comparison, the calculated 
binding energy of 4 He, 27.408 MeV, obtained by using the same no-core approach in the 
8hQ space |§, was added to the diagonal matrix elements of our hamiltonian. Note that 
by diagonalizing the hamiltonian in Table DL we get the same excitation energies as those 
from the 6hQ calculation given in the third column of Table IV of Ref. ||. We also present 
the effective hamiltonian obtained in the same way from a 4hQ calculation. Note that the 
dimension of the 4hfl calculation is more than an order of magnitude smaller than that of 
the 6hQ calculation. The effective interaction used in this calculation was obtained from 
the 6hQ multi-valued interaction by leaving out the set corresponding to the 6hfl space. 
Let us mention that the change in the low-lying eigenenergies in the two calculations is not 
substantial and the ordering of levels is identical. We observe that our calculated matrix 
elements differ in some cases from the phenomenological ones. Let us point out, however, 
that our matrix elements provide a better description of 6 Li states than those of Cohen- 
Kurath. This can be understood as the latter matrix elements were fitted to a large number 
of nuclei across the entire Op-shell. 

Our primary aim is to derive Op-shell effective electromagnetic operators. In the full-space 
calculation we employed the one-body E2 and Ml operators 
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»=i 

+e n E(|-^K 2 F( 2 )(9 i ), (5a) 



i=l 

A 



t{m1) = vJ^e[(I+^)w+^ 



+(£-*«)(</fl< + , (5b) 

with the free nucleoli charges e p = e, e n = and free nucleon (/-factors gf = 1, gf = 0, 
gv = 5.586, and g* = —3.826. The P-space operators are constructed by the application 
of Eq. (|4]). We calculate the P-space T( E2 ) operator and separately the orbital and spin 
parts of the Ml operator. The calculation is performed in the m-scheme and subsequently 
transformed to the J, T basis. To get all the reduced matrix elements, full-space calculations 
with Mj = and Mj = 1 must be done. In Table || we present selected matrix elements of 
pieces of (j5j), namely the operators 

E2 = VlFEG + t^rfY^iQi) , (6a) 

i=l 

A 

M\ = J2(l + t Z i)k, (6b) 
i=i 

A 

Ms = + t zi ) Si . (6c) 
i=i 

These matrix elements are reduced in J and for T z = 0. The second column shows the 
matrix elements of the effective operators, as obtained from Eq.(^j) and the procedure out- 
lined above. Note that these operators, when used with the eigenvectors of the effective 
hamiltonian obtained from (^), whose matrix elements are shown in Table |, give the same 
mean values and transition rates as the original one-body operators (@), when used with the 
OhQ dominated eigenvectors of the 6hQ calculation. Also note that the effective operators 
are two-body operators unlike the full-space original operators. 

Let us first discuss the E2 operator. In the third column of Table [TT] the reduced matrix 
elements of the operators (|6|), evaluated in the P-space, are shown for comparison. We 
observe, that there is a striking difference in the renormalization of the isoscalar and isovector 
matrix elements of (E2) e fj. The former are much larger in magnitude than the latter in 



comparison with the unrenormalized values of the operator (|6a|) . Apparently, there is no 



chance to approximate the effective operator as ( |6a| ) multiplied by some effective charge. 
Instead, it is possible to mimic the mentioned isoscalar-isovector effect by approximating 
the effective operator as a combination of one-body proton and neutron operators with 
different effective charges, e.g., 



(E2) eff ^eP ffV ^E(| + ^)^( 2 )(e i ) 

8=1 

+< ffv /I FE(i-^K 2 ^ 2) (0 l ), (7) 



i=l 
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where only valence nucleons contribute in the sums. A better approximation may be 
obtained when the effective charges become j- dependent, e.g., e cSij(i\0\j)a\aj in the 
second-quantization form. We calculated the effective charges from the reduced matrix ele- 
ments of the appropriate operators by a least-square fit. The resulting j-dependent, as well 
as j-independent, effective charges are presented in Table [ITT], and the corresponding reduced 
matrix elements are shown in the fourth and fifth columns of Table ||. First, we observe 
that this kind of approximation works very well. Moreover, the pure two-body matrix ele- 
ments, which cannot be reproduced by an approximation of the type (0), are almost a factor 
of ten smaller than the largest one-body matrix elements. Second, the calculated effective 
charges e^ s = 1.527e, e° ff = 0.364e are close to the phenomenological ones mentioned in the 
introduction. Third, the j-dependence of the effective charges is rather moderate. 

From the phenomenological studies it is well known that the magnetic dipole transitions 
and moments can be in most cases described, at least in light nuclei, by using the operator 
( pE|) with small modification of the g-factors. Also the effective orbital and spin operators 
obtained in our study are much less renormalized, when compared to the starting one-body 
operators, than the electric quadrupole operator. It is reflected in the second and third 
parts of Table |I|. The isoscalar-isovector effect is much smaller, and, in the case of the 
spin operator it is almost non-existent. Unlike the case of the quadrupole operator, the 
effective dipole operator matrix elements are, on average, smaller in comparison with one- 
body ones. A perhaps surprising result is, however, the observation that the orbital part is 
more renormalized than the spin part and, moreover, a neutron orbital part is generated with 
an effective ^-factor gf e ^ = 0.085. The proton ^-factor is about 10% quenched gf e g = 0.907. 
The spin part is quenched by about 6%, e.g., g^ cS = 0.937#p + 0.001g°, and from the isospin 
symmetry = 0.937g" + O.OOlgP. As in the quadrupole operator case, the j-dependence 
is also moderate for the magnetic dipole operator effective charges, but the pure two-body 
matrix elements are relatively smaller. 

The j-independent effective charges extracted in the same way from the 4hQ calculation 
are also presented in Table |T|. We observe that the renormalization is smaller in this 
case. This reflects the fact that, for example, the E2 transition rates obtained in the 4hQ 
calculation are weaker than those calculated in the 6hQ space. Our observation here is that 
the effective hamiltonian converges more rapidly than the E2 operator with respect to the 
full-space size change. 

To quantify the two-body content of the effective Op-shell electromagnetic operators, we 



evaluate the quantity R = yJ2i<j(O e s,ij — o h- c s,ij) 2 / J2i<j(O e s,ij) 2 , where the Oy denote 
the matrix elements between the Op-shell two-body states i and j. In this way we estimate 
the part of the effective operators, which cannot be expressed as a combination of one-body 
operators. When using one-body operators with j-dependent charges, we obtain a two-body 
contribution of 10.1% for E2, 3.8% for Ml, and 3.2% for Ms, respectively. For the one-body 
operator with j-independent effective charges the two-body contributions are 12.3% for E2, 
5.3% for Ml, and 4.4% for Ms, respectively. Clearly, the magnetic dipole operators are better 
approximated by combinations of one-body operators. 

In conclusion, we have shown that model-space truncation is sufficient to generate op- 
erator renormalization, which is characterized by effective charges compatible with those 
used in the phenomenological applications. We have found that the isoscalar and isovector 
parts of the operators are renormalized differently, particularly in the case of the electric 



5 



quadrupole operator. This difference in renormalization is the source of a non-zero neutron 
effective charge. These findings are based on a no-core 6hQ calculation for 6 Li with a multi- 
valued starting-energy-independent two-body interaction derived microscopically from the 
Reid 93 nucleon-nucleon potential, from which the Ohil Op-shell effective hamiltonian, elec- 
tric quadrupole, and magnetic dipole operators were constructed. The obtained effective 
operators are two-body operators. We have shown, however, that they may be well approx- 
imated by one-body operators. Their two-body content is about 10% for E2 and about 4% 
for Ml opertators, respectively. We also studied the dependence of the renormalization on 
the size of the full space. We observed a non-negligible difference between the effective E2 
charges extracted from the 6hQ calculation in comparison to those obtained from the AMI 
calculation. On the other hand, the changes in the effective hamiltonians obtained in the 
two calculations are less pronounced. 
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TABLES 



(2ji2j 2> JT 


H 


2j 3 2j 4 , JT) 


Eff 


Eff-4 


CK 6-16 


(11,01 


H 


11,01) 


6.772 


7.165 


4.88 


(11,01 


H 


33,01) 


-2.756 


-3.201 


-5.32 


(33,01 


H 


33,01) 


0.493 


-0.970 


0.52 


(11,10 


H 


11,10) 


3.999 


4.202 


0.28 


(11,10 


H 


13,10) 


-0.776 


-0.988 


-1.39 


(13,10 


H 


13,10) 


-0.788 


-1.333 


-2.64 


(11,10 


H 


33, 10) 


2.086 


2.464 


1.09 


(13,10 


H 


33, 10) 


-4.107 


-4.708 


-4.02 


(33, 10 


H 


33, 10) 


0.815 


-0.707 


0.12 


(13,11 


H 


13, 11) 


5.780 


5.308 


4.76 


(13,20 


H 


13,20) 


0.199 


-1.366 


-0.32 


(13,21 


H 


13,21) 


4.303 


3.514 


2.76 


(13,21 


H 


33,21) 


1.377 


1.613 


2.21 


(33,21 


H 


33,21) 


2.694 


1.162 


2.61 


(33,30 


H 


33, 30) 


-1.842 


-3.771 


-3.42 


TABLE I. 


The Op-shell effective 


hamiltonian 


matrix elements, in MeV, 


obtained from 6hQ, 



second column, and from AHQ, third column, calculation for 6 Li. The calculated binding energy 
of 4 He, 27.408 MeV, obtained by using the same method in the 8h£l space ||, was added to 
the diagonal matrix elements in order to make a meaningful comparison with the Cohen-Kurath 



phenomenological matrix elements [14], presented in the fourth column 



S 



(O) 


Eff 


ob 


ob eff-j 


ob eff 


(11, 10|E2|13, 20) 


5.916 


2.739 


5.539 


5.179 


(11,10|E2|13,21) 


-3.101 


-2.739 


-3.256 


-3.184 


(13, 20 E2 13,21) 


2.137 


2.092 


2.321 


2.432 


(13, 21[E2|13, 21) 


-3.612 


-2.092 


-3.605 


-3.956 


(13, 10|E2|33,21) 


-2.481 


-1.937 


-2.303 


-2.251 


(13, 10 E2 33,30} 


-6.295 


-3.240 


-6.554 


-6.128 


(33,01 E2|33, 21} 


-3.942 


-2.236 


-3.854 


-4.229 


(33, 21 E2 33, 30} 


4.122 


3.742 


4.152 


4.350 


(11, 10 E2 11,10} 


0.501 


0.0 


0.0 


0.0 


(11,10|E2|33,30) 


0.649 


0.0 


0.0 


0.0 


(11, 01|M1|11, 10) 


-0.902 


-1.155 


-0.912 


-0.949 


(11, 10|M1|11, 10) 


1.560 


1.633 


1.585 


1.621 


(13, 11 Ml 33, 10) 


-0.439 


-0.646 


-0.509 


-0.531 


(13,21 Ml 33,21) 


-0.603 


-0.646 


-0.586 


-0.641 


(33,01 Ml|33, 10) 


-1.098 


-1.291 


-1.096 


-1.061 


(33, 30 Ml 33, 30) 


3.122 


3.055 


3.106 


3.032 


(11,01 Ml 33,10) 


-0.100 


0.0 


0.0 


0.0 


(ll,01|Ms|ll,10) 


0.289 


0.289 


0.280 


0.270 


(11, 10|Ms|ll, 10) 


-0.336 


-0.408 


-0.360 


-0.383 


(11, 10 Ms 13,20} 


0.804 


0.913 


0.829 


0.856 


(ll,10Ms|13,21) 


-0.819 


-0.913 


-0.839 


-0.855 


(33,01 Ms 33, 10} 


-0.618 


-0.646 


-0.620 


-0.604 


(33, 30 Ms 33, 30} 


1.460 


1.528 


1.476 


1.433 


(11,01 Ms 33,10} 


0.041 


0.0 


0.0 


0.0 



TABLE II. Selected reduced matrix elements of the proton E2, in h/mQ,, and Ml, in /^n, 
operators. Here, (O) = (2j'i2j2, J\Ti 1 2j3 2 j'4 , J3T3}. In the second column the Op-shell effective 
operator matrix elements, obtained from the 6hU calculation for 6 Li, are presented. The third 
column shows the corresponding proton one-body (ob) operator matrix elements. The fourth 
and fifth columns display the matrix elements of the one-body operators with j-dependent and 
j-independent effective charges, respectively. These operators are combinations of one-body proton 
and neutron operators. 
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e p 13 

X o 
2 2 


2 2 


p p 

3 3 
2 2 


033 
22 


e 1 ?! 

1 J 

2 2 


e 1 1 
2 2 


iLil 


1.606 


0.417 


1.417 


0.307 






Ml 


0.848 


0.060 


0.933 


0.084 


0.880 


0.090 


Ms 


0.914 


-0.006 


0.963 


0.003 


0.925 


-0.043 




p p 


e eff 


p p 


p n 






E2 


1.527 


0.364 


1.302 


0.244 






Ml 


0.907 


0.085 


0.931 


0.063 






Ms 


0.937 


0.001 


0.953 


-0.003 







TABLE III. Effective charges of the proton quaelrupole, magnetic orbital and magnetic spin 
operators, derived by least-square fits to the corresponding Op-shell effective operators obtained 
from the QhCl calculation for 6 Li. Both the j-dependent and j-independent effective charges are 
shown. Also, the j-independent effective charges obtained in the same way from the AfiQ, calculation 
are presented in the last two columns. For the definition of the effective charges see Eq. (|7|). 
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